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SUMMARY 


Modern  numerical  methods  for  the  optimization  of  large  dis¬ 
cretized  systems  are  now  we  1 1  developed  and  highly  efficient  in 
the  case  of  thin  \v’ ailed  elastic  structures  mcaeled  by  finite 
elements.  However,  this  is  not  yet  true  for  structures  whose 
components  are  subject  simultaneously  to  bending  and  extension 
loads.  In  this  paper,  the  idea  of  Generalized  Optimality  Cri¬ 
terion  (GOC),  set  forth  in  previous  final  scientific  reports 
for  bar,  membrane  and  pure  bending  elements,  is  extended  to  deal 
with  general  beam  and  flat  shell  elements.  The  modifications 
brought  to  the  GOC  result  in  explicit  approximations  for  the 
behavior  constraints  that  are  still  correct  up  to  the  first  order, 
but  that  exhibit  a  more  comp,  lex  algebraic  form.  Indeed  these 
explicit  expressions  are  no  longer  merely  linear  in  the  reciprocal 
design  variables.  However  they  continue  to  be  addi lively  sepa¬ 
rable  and  therefore,  dual  methods  remain  fully  applicable,  just 
as  in  the  original  statement  of  the  GOC  approach.  Numerical 
examples  will  be  offered  to  demonstrate  the  efficiency  of  the 
method  presented. 
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INTRODUCTION 


The  optimum  design  of  any  significant  structure  is  the  result 
of  a  delicate  compromise  between  many  complex  factors.  Some  art- 
rational  and  can  be  quantified,  such  as  the  strength  of  t fit- 
structure,  its  natural  frequencies,  its  weight,  ...  Some  are 
just  as  rational  but  are  difficult  to  quantify,  such  as  the  ex¬ 
perience  of  the  designer  in  a  given  technology.  Some  others  are 
much  less  rational,  like  styling,  but  are  just  as  important  for 
the  final  goal  of  the  process,  which  is  the  marketing.  Naturally 
a  good  designer  considers  structural  optimization  as  a  technique 
that  should  take  into  account  all  possible  aspects  of  the  design. 
Consequently  the  designers  are  often  reluctant  to  the  concepts 
of  structural  optimization  developed  in  connection  with  finite 
element  programs. 

However  a  more  detailed  examination  of  the  design  process 
leads  to  isolate  a  phase  that  appears  frequently,  during  which 
the  shape  of  the  structure  is  more  or  less  frozen  and  the  problem 
is  limited  to  giving  adequate  dimensions  to  the  various  members. 
Such  a  situation  is  often  encountered  in  the  aerospace,  naval 
or  automobile  industries,  where  the  external  shapes  are,  to  a 
large  extent,  dictated  by  aero-  or  hydrodynamic  considerations, 
or  by  styling,  while  internal  forms  are  often  determinated  by 
various  other  non  structural  requirements.  If  the  ultimate  goal 
of  the  designer  can  be  identified  as  corresponding  to  the  minimi¬ 
zation  of  an  explicit  function  of  the  member  sizes,  and  if  the 
limitations  on  the  design  can  be  defined  as,  eventually  implicit, 
functions  of  the  member  sizes  too,  such  as  displacements,  stresses, 
eigenfrequencies,  etc...,  then  the  problem  is  tractable  by  auto¬ 
matic  algorithms.  They  allow  the  designer  to  speed  up  signifi¬ 
cantly  this  part  of  the  design  process  and  to  explore  more  sys¬ 
tematically  the  various  feasible  designs. 

The  optimization  problems,  which  in  fact  should  be  called 
automatic  sizing  problems,  are  especially  crucial  when  complex 
structural  forms  are  involved,  when  flexural  forces  cannot  be 


2  . 

neglected,  and  when  composite  materials  such  as  reinforced 
resins  are  employed.  In  these  cases  it  becomes  difficult,  if 
not  impossible,  for  the  designer  to  have  an  intuitive  under¬ 
standing  of  the  structural  mechanics  that  is  sufficient  to  1 c  id 
to  op t i mal  s  i  z i ng  of  the  various  members.  Furthermore,  the  de¬ 
signer  is  most  of  the  time  unable  to  take  into  account  global 
constraints  in  the  structure,  like  global  flexibility,  restric¬ 
tion  on  displacements,  frequencies  of  vibration,  global  buckling 
modes,  etc...  It  is  only  possible  to  verify  a  posteriori  that 
such  constraints  are  satisfied.  Again  these  global  constraints 
become  more  important  in  the  context  of  highly,  indeterminate 
structures.  In  the  aerospace  industry,  the  necessity  of  desi¬ 
gning  high  performance  structures  has  motivated  significant  re¬ 
search  efforts  to  derive  algorithms  permitting  a  rapid  and  sys¬ 
tematic  exploration  of  the  design  space  to  determine  the  optimum 
material  utilization. 

It  is  worth  pointing  out  that  optimization  methods  should 
be  considered  as  especially  useful  in  the  preliminary  design 
phase.  Using  them  when  the  design  is  practically  frozen,  with 
the  hope  of  an  ultimate  improvement,  is  often  disappointing. 

This  i  due  to  the  fact  that  the  optimization  of  a  detailed  design 
implies  the  formulation  of  a  large  number  of  constraints,  some 
of  which  are  not  easily  quantified.  At  the  preliminary  design 
stage,  however,  the  constraints  are  usually  more  global  and 
therefore  more  easily  handled  by  the  available  formulations. 

The  structural  optimization  problem  considered  in  this  report 
consists  of  the  weight  minimization  of  a  finite  element  model 
with  fixed  geometry  and  material  properties.  The  design  variables 
are  taken  as  the  transverse  sizes  of  the  structural  members, 
namely,  the  cross-sectional  areas  of  bar  and  beam  elements  and 
the  thicknesses  of  membrane,  plate  and  flat  shell  elements. 

The  mathematical  programming  problem  to  be  solved  has  the  fol¬ 
lowing  form  : 


minimi ze 

n 

W  *  .  I  .  1.  a  . 

i*l  i  i 

(  1  ) 
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The  a^'s  denote  the  n  design  variables  .  They  correspond  to 
member  sizes  of  either  individual  finite  elements,  or,  if 
design  variable  linking  is  used,  of  groups  of  finite  elements. 
The  structural  weight  W  is  a  linear  objective  function,  because 
the  £  ^  '  s  are  constant  coefficients  representing  the  specific 
weight  times  either  the  element  length  (bars  and  beams)  or  the 
element  area  (membranes,  shearpanels,  plates  and  flat  shells). 
The  inequalities  (2)  are  the  behavior  constraints,  which  impose 
limitations  on  quantities  describing  the  structural  response, 
for  example,  the  stresses  and  displacements  under  multiple 
loading  cases,  the  natural  frequencies,  the  buckling  loads, 
etc...  The  design  variables  are  also  subjected  to  the  side 
constraints  (3),  where  a^  and  a^  are  lower  and  upper  limits  that 
reflect  fabrication  and  analysis  validity  considerations. 

Standard  minimization  techniques  can  be  applied  to  the  non¬ 
linear  programming  problem  (1-3).  However  this  problem  exhibits 
some  characteristics  that  make  it  complicated  when  practical 
structural  design  applications  are  considered.  The  essential 
difficulty  arises  from  the  implicit  nature  of  the  behavior 
constraints  (2),  in  that  their  precise  numerical  evaluation  for 
each  particular  design  requires  a  complete  finite  element  ana¬ 
lysis.  Since  the  solution  scheme  is  iterative,  it  involves  a 
large  number  of  structural  reanalyses.  Therefore  the  computa¬ 
tional  cost  often  becomes  prohibitive  when  large  structural 
systems  are  dealt  with.  However  a  powerful  design  procedure 
has  now  emerged,  which  consists  in  replacing  the  initial  problem 
with  a  sequence  of  simple  explicit  problems.  In  the  next  sec¬ 
tion  this  approach  will  be  briefly  reviewed  by  restricting  the 
formulation  to  thin  walled  structures  idealized  by  bar/membrane 
elements.  A  much  more  detailed  presentation  can  bo  found  in  a 
previous  report  [  1  ] .  It  will  be  shown  that  the  behavior  cons¬ 
traints  can  be  approximated  either  by  using  virtual  load  consi¬ 
derations  (optimality  criteria  approach)  or  by  using  first  order 
Taylor  series  expansion  with  respect  to  the  reciprocal  design 
variables  (mathematical  programming  approach).  Applying  a  dual 
solution  scheme  to  each  explicit  problem  generated  in  sequence 
naturally  introduces  the  concept  of  a  generalized  optimality 


criterion. 


S ub s eq ue n t 1 y s e e t i on s  3  and  4  wilL  be  concerned  with  struc¬ 
tural  models  that  are  capable  of  carrying  flexural  forces.  For 
beams  and  plates  in  pure  bending,  adequate  intermediate  varia¬ 
bles  can  be  selected,  in  terms  of  which  high  quality  explicit 
approximations  for  the  behavior  constraints  can  still  be  gene¬ 
rated  by  linearization.  The  idea  of  generalized  optimality 
criterion  remains  fully  valid  and  it  keeps  its  interpretation 
in  terms  of  energy  densities  in  the  structural  members.  Section 
3  is  a  summary  of  results  presented  in  a  previous  report  (ref. 

[2  j,  section  6). 

The  essential  problem  is  addressed  in  section  4.  It  consists 
of  the  establishment  of  the  generalized  optimality  criterion 
approach  in  the  general  case  where  the  structural  members  work 
both  in  extension  and  flexion  (beam  and  flat  shell  elements). 

For  displacement  constraints,  it  is  no  longer  possible  to  select 
a  suitable  intermediate  variable  for  the  linearization  process, 
however,  the  virtual  load  procedure  permits  obtention  of  high 
quality,  first  order  explicit  approximations  of  the  behavior 
constraints.  It  will  also  be  shown  how  to  proceed  for  stress, 
frequency  and  buckling  constraints. 

The  explicit  approximations  still  exhibit  a  separable  alge¬ 
braic  form  and  therefore  dual  methods  remain  applicable.  However 
some  difficulties  might  happen  due  to  the  lack  of  convexity  in 
the  explicit  subproblera.  Section  3  will  describe  how  to  solve 
the  explicit  subproblems  by  using  dual  methods.  Finally  seine 
applications  will  be  offered  in  section  6  to  illustrate  the  power 
and  generality  of  the  approach  presented. 


2  .  GENERALIZED  OPTIMALITY  CRITERION  FOR  THIN  WALLED  STRUCTURES 

This  section  summarizes  some  results  obtained  in  a  previous 
work  (  1  J  for  structural  models  made  up  of  bar,  membrane  and 
shear  panel  elements,  which  are  quite  adequate  for  idealizing 
thin  walled  structures  subjected  mainly  to  extension  loading. 
For  this  class  of  finite  element  models,  the  structural  stif¬ 
fness  matrix  exhibits  a  linear  form  in  the  design  variables  : 


K  = 


K  . 
1 


i=  I 


a  . 
i 


( '•  ) 


where  K^,  a  matrix  of  constants,  represents  the  stiffness  matri 


of  the  it  It  element  when  a  ^  =  ( 


For  simplicity,  the  following 
discussion  is  restricted  to  problems  involving  constraints  on 
static  stresses  and  displacements,  in  which  case  the  behavior 
constraints  (2)  can  be  written 


Va) 


u  . 
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u . (a)  >  0 
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(  '>  ) 


where  u.  denotes  an  upper  bound  to  a  response  quantity  Ujf.i) 
(stress,  nodal  displacement,  relative  displacement). 


Most  of  the  optimality  criteria  approaches  (e.g.  I  3  1  ),  a  fa- 
well  as  the  generalized  optimality  criterion  (GOC)  set  forth  in 
Ret.  |  I  ] ,  use  the  virtual  load  technique  to  generate  explicit 
approximations  of  the  stress  and  displacement  constraints.  In¬ 
troducing  a  virtual  load  vector  conjugated  to  the  response  quan¬ 
tity  Uj  (unit  load  for  a  nodal  displacement),  it  follows  that 
u j  can  be  expressed  as  the  sum  of  the  contributions  of  each  fi¬ 
nite  element  : 


T  n  Ci  j 

uj  =  q  K  <ij  =  ;Ii  — 1 


i=l  a . 
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with 


T  -  2 

c..  =  (q  K.  q.)  a. 

l  j  l  i 
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Ln  these  expressions  q  and  are  respectively  the  real  and 

virtual  displacement  vectors  and  K  .  is  the  element  s  t  i  !  f  n e  , 

matrix  appear:  hr  in  (■'*).  It  can  be  seen  from  (7)  that  tlu-  •  i>e!- 

ficients  e.  .  are  related  to  the  virtual  strain  eiu-rvv  densities 
i  J 

in  the  structural  members.  The  c..'s  are  constant  coefficients 

i  J 

in  the  case  of  a  statically  determinate  structure,  so  that  ((>) 
represents  then  the  exact  explicit  form  of  the  response  quan¬ 
tity  Uj.  In  the  case  of  a  statically  indeterminate  structure, 
the  c^j's  depend  implicitly  on  the  design  Vdtiabler,  because 
structural  redundancy  produces  redistribution  of  the  internal 
forces  when  the  member  sizes  are  modified.  Therefore  the  loi- 
1  owing  explicit  constraints 


•\j  -  n  i  j 

h  .  ( a )  u  .  -  .  >i  .  - — J-  >  0  18) 

J  J  i= 1  a. 

constitute  in  general  approximate  forms  of  the  original  cons¬ 
traints  (5).  As  shown  in  Ref.  j  1  ] ,  the  basic  idea  in  the  op¬ 
timality  criteria  approach  can  be  viewed  as  transforming  the 
initial  implicit  problem  into  a  sequence  of  explicit  subproblem 
Each  explicit  problem  results  from  replacing  the  behavior  cons¬ 
traints  (2)  by  their  approximate  forms  (8). 


On  the  other  hand,  the  mathematical  programming  approach  to 
structural  optimization,  after  a  period  of  unefficiency,  has 
finally  evolved  into  a  powerful  and  now  well  established  design 
procedure  which  is  also  based  upon  explicit  approximations  o I 
the  behavior  constraints  [4,  3,  6  ).  The  key  idea  is  to  linea¬ 
rize  the  behavior  constraints  witli  respect  to  the  reciprocal 
design  variables 


Justification  for  this  change  of  variables  lies  in  the  fact 
that  the  constraint  surfaces  can  be  shown  to  be  very  shallow 
and  close  to  planes  in  the  reciprocal  design  variable  space. 
Therefore  the  linearized  forms  of  the  constraints  are  usually 
high  quality  approximations.  They  are  obtained  by  using  first 
order  Taylor  series  expansion  in  terms  of  the  reciprocal  va¬ 
riables  x.  : 
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lij  (x) 


l  u 


U .  -  luT  *  .  >: ,  ( 

J  .]  1  =  1 
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where  the  u|>  pcrscript  denotes  quantities  evaluated  it  tin- 

actual  design  point  x°,  where  the  structural  analysis  is  j  > .  •  r 
formed.  Note  that  the  finite  element  analysis  capability  •mi:,  t 
include  auxiliary  sensitivity  analyses  for  evaluating  the  lirsi 
partial  derivatives  of  the  response  quantities.  Most  often  the 
well  known  pseudo-loads  technique  is  employed  |  7  ]  . 

It  has  been  shown  in  a  previous  report  |  I  |  that  the  explicit 
approximations  of  the  behavior  constraints  used  in  both  the 
optimality  criteria  and  mathematical  programming  approaches  ( 1.  q  t, 

8  and  10,  respectively)  are  identical.  Indeed  the  virtual  si  rain 
energy  densities  employed  in  the  optimality  criteria  appro¬ 

aches  are  nothing  else  that  the  gradients  of  the  response  quan¬ 
tities  with  respect  to  the  reciprocal  variables  : 


c  .  . 
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i  J  3  x  . 
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Furthermore  the  definition  of  the  c..'s  following  from  the  vir- 

i  J 

t  u  a  1  load  technique  (see  E  q  .  7)  clearly  indicates  that 
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Therefore  (10)  can  be  rewritten 


h  j  (  x  )  u  j 
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(it) 


which  is  equivalent  to  (8)  when  recast  in  terms  of  the  direct 
variables  a^.  It  is  thus  apparent  that  a  unified  approach  to 
structural  weight  minimization  of  finite  element  systems  has 
emerged,  which  consists  in  replacing  the  initial  problem  (1-3) 
with  a  sequence  of  explicit  approximate  -  or  linearized  -  pro¬ 
blems  of  the  following  form  : 


minimize 
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The  GOC  statement  results  from  writing  the  KL'HN-TUCKEK  optimal  i  tv 
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Of  course  the  lagrangian  multipliers  must  be  uonne native ,  more 
precisely,  they  must  satisfy  the  complementary  conditions  : 
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(20- 

21), 

an  interesting 

approach  is  to 

resort 

to  dual  methods, 
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dual 

problem 

has  been  found,  the  corresponding 

optimal  des 

i  n 
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relations  (17- 

19).  Note  that 

the  design 

variables  can  be 

separated  into  a  group  of  active  (or  free)  variables  (see  l.q  .  17) 

and  a  group  of  passive  (or  fixed)  variables  (see  Eqs.  18  and  19). 
This  subdivision  into  active  and  passive  design  variable  groups 
is  classical  in  the  optimality  criteria  approaches  |  3,  9  ]  .  It 

corresponds  to  the  fact  that  the  dual  space  -  i.e.  the  space  of 
the  lagrangian  multipliers  r^  -  is  partitioned  into  several  re¬ 
gions  separated  by  planes  across  which  the  second  derivatives 
of  the  dual  function  are  discontinuous  [6,  8  ]  . 


The  whole  process  of  combining  the  linearization  of  the 
behavior  constraints  with  respect  to  the  reciprocal  design 
variables  and  a  dual  solution  scheme  can  be  viewed  as  a  gene¬ 
ralization  of  the  optimality  criteria  approaches.  It  is  impor¬ 
tant  to  mention  that  this  basic  approach  of  converting,  the 
initial  problem  into  a  sequence  of  explicit  subproblems  is  now 
widely  recognized  [6,  10  ]  and  it  is  routinely  employed  for 

large  scale  industrial  applications  [  11  ]  . 


3  .  P  U  R  L  B  V.  N  D  [  NG  F  L  K  M  ENTS 

In  this  section  attention  is  focused  on  discretized  models 
made  up  of  pure  beam  and  plate  elements  subjected  to  flexural 
loads  only  (for  more  details,  see  section  6  of  Ref.  |2  J). 

The  stiffness  matrix  of  such  a  bending  element  is  usually  not 
merely  proportional  to  its  cross-sectional  size  and  therefore 
the  optimization  strategy  reviewed  in  the  previous  section  must 
be  modified.  The  way  to  deal  with  a  beam  element  subjected  to 
uniaxial  bending  depends  upon  the  relationship  between  the  prin¬ 
cipal  moment  of  inertia  I  and  the  cross-sectional  area  a.  A 
wide  variety  of  situations  is  taken  into  consideration  by  adop¬ 
ting  the  following  relation  : 

I  =  c  aP  (22) 


where  c  is  a  constant  that  depends  only  on  the  shape  of  the 
beam  cross-section  and  p  is  a  positive  number. 


Most  of  the  time  p  is  taken  as  -n  integer  number,  equal  to 
1,  2  or  3.  The  case  p  =  1  corresponds  to  thin  walled  beams, 
for  example,  sandwich  beams,  pipes  with  fixed  diameter  and  va¬ 
riable  thickness,  etc...  The  GOC  approach  of  section  2  remains 
then  fully  applicable,  since  the  stiffness  matrix  continues  to 
be  linear  in  the  design  variables.  The  case  p  =  2  is  that  of 
beams  with  uniformly  varying  cross-section.  The  shape  of  the 
cross-section  is  kept  unvariant  while  its  area  is  modified  during 
redesign  (dilatation  or  contraction).  Finally  the  case  p  =  3 
is  concerned  with  beams  having  full  cross-section  whose  height 
varies  while  other  sizes  are  fixed.  For  a  beam  subjected  to 
pure  bending,  the  flexural  rigidity  is  proportional  to  the  mo¬ 
ment  of  inertia  and  therefore,  in  a  finite  element  context,  the 
structural  stillness  matrix  exhibits  the  following,  explicit 
form  in  terms  of  the  cross-sectional  areas  : 
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1  =  1  i 


a?  K. 


where  each  matrix  is  independent  of  the  design  variables 
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With  regard  to  plate  elements  subjected  to  pure  bending, 
two  cases  must  be  distinguished.  The  first  case  is  that  of 
sandwich  plates  with  constant  core  thickness.  The  sheet  thick¬ 
nesses  constitute  then  the  design  variables.  Consequently  the 
stiffness  matrix  continues  to  depend  linearly  on  the  design 
variables  and  the  GOC  relations  (17-21)  remain  fully  applicable 
The  second  case  is  concerned  with  full  plates  with  variable 
thickness.  The  stiffness  is  then  proportional  to  the  cube  ol 
the  thickness  and,  in  an  assembling  of  plate  elements,  the  re¬ 
lation  (22)  must  be  chosen  with  p  =  3. 

By  assuming  that  the  structural  discretization  is  made  up 
entirely  of  elements  of  the  same  type,  the  stiffness  matrix 
exhibits  the  form  (23),  where  p  takes  on  the  same  value  fru 
each  member.  In  these  circumstances,  the  GOC  can  be  derived 
just  as  in  the  case  of  thin  walled  structures,  by  adopting  a 
change  of  variables  tending  to  reduce  the  nonlinear  character 
of  the  constraints  : 


a . 

l 

The  next  step  is  to  linearize  the  constraints  with  respect  to 
the  new  variables  x^,  which  requires  gradient  evaluation  (see 
E  q  .  10).  Restricting  again  the  discussion  to  stress  and  dis¬ 
placement  constraints,  it  is  easily  shown  (see  Ref.  I  2  |)  that 
the  first  order  Taylor  series  expansion  (10)  reduces  to  the 
form  (13),  or,  when  written  in  terms  of  the  direct  variables  : 


h  .  ( a )  r  u  .  -  .  . 

J  J  1=1 
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(23) 


The  c^j  coefficients  can  be  interpreted  as  the  gradients  of  the 

response  quantities  with  respect  to  the  intermediate  variables 

x^  defined  in  (24),  but  they  can  also  be  related  to  the  virtual 

strain  energies  e. .  in  the  structural  members  : 

i  J 
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i  j 
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where,  by  definition. 


e  .  .  =  q  *  K  -  o  .  (  .’  7  ) 

>  1  i  l 

In  this  connection  it  should  he  recognized  that  virtual  load 
considerations  could  directly  be  employed  to  derive  the  explicit 
approximations  (25),  instead  of  resorting  to  first  order  Taylor 
series  (see  next  section). 


The  optimality  criteria  equations  are  very  similar  to  (17- 
19),  provided  that  care  is  taken  of  the  exponent  p  appearing, 
in  (25).  The  basic  redesign  relations  for  the  active  desi gn 
variables  must  read  as  follows  : 
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(2b) 


where  it  is  understood  that  the  dual  variables  r.  (i.e.  the 
lagrangian  multipliers)  must  satisfy  the  complementarity  con¬ 
ditions  (see  Eqs.  20  and  21).  A  physical  interpretation  of  the 
optimality  criterion  is  obtained  by  introducing  the  virtual 
strain  energy  densities  per  unit  weight  : 
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where  e..  is  defined  in  (27).  in  terms  of  the  <  .  .  ’  s  ,  the  o  p  t  i 
tj  i) 

mality  criterion  (28)  takes  the  "energetic"  form 
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In  the  special  case  where  only  one  displacement  constraint 
is  specified,  the  optimality  criterion  states  that  the  virtual 
strain  energy  density  must  be  the  same  in  each  element.  In 
this  simple  case,  it  is  possible  to  solve  analytically  the  ex¬ 
plicit  problem  and  to  derive  explicit  redesign  relations  in  terms 
of  known  quantities.  The  active  design  variables  can  be  shown 
to  be  given  by 
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while  the  remaining  passive  variables  are  fixed  to  an  upper  or 
a  lower  limit  (  uq  denotes  the  contribution  of  these  passive 
variables  to  the  displacement  constraint  u  <  u).  It  is  worth 
mentioning  that  (31)  is  well  suited  for  the  design  o  f  plates 
in  bending  with  a  single  displacement  constraint.  Since  then 
p  =  3  the  redesign  relation  (31)  involves  the  fourth  root  of  the 
coefficients  c  ^ ,  rather  than  the  third  root  as  employed  in 
Ref.  [12  1  on  an  intuitive  basis.  Note  also  that  by  taking 
p  =  1  in  (31),  conventional  redesign  relations  are  recovered, 

which  were  devised  for  trusses  [3  ]  ,  sandwich  beams  (  I  i  j  .etc... 

It  can  be  concluded  that  the  generalized  optimality  criteria 
approach  can  easily  be  extended  to  deal  with  pure  bending  ele¬ 
ments  by  defining  adequate  intermediate  variables.  High  quality 
explicit  approximations  of  the  behavior  constraints  can  s  t i  1  1  be 
generated  and  the  resulting  GOC  keeps  its  interpretation  in  terms 
of  energy  densities  ir  the  structural  members.  Seeking  the  de¬ 
sign  variables  that  satisfy  the  GOC  at  each  redesign  stage  can 
still  be  achieved  efficiently  by  resorting  to  dual  methods,  be¬ 
cause  the  explicit  approximate  problem  remains  separable  and 
strictly  convex  when  expressed  in  the  intermediate  design  vari¬ 
ables  x  ^  . 


4.  FLEX  l  ON-K  XTEN  S ION  ELEMENTS 


When  flexion  and  extension  loadings  act  simultaneously  witli 
comparable  intensity  at  the  element  level,  the  definition  (21) 
of  the  stiffness  matrix  can  no  longer  characterize  the  structu¬ 
ral  model  with  sufficient  accuracy.  To  help  fix  ideas,  consi¬ 
der  a  flat  shell  element  made  up  of  a  membrane  and  a  plate  sta¬ 
cked  together.  The  stiffness  matrix  of  such  a  Mat  shell  ele¬ 
ment  exhibits  the  form 


K.  . 
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♦  a3  K.(3) 
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(  32  ) 


where  and  xf3^  are  constant  matrices.  As  a  result,  in 

t  l 

the  GOC  approach,  if  the  constraints  are  linearized  with  respect 
to  the  reciprocal  design  variables  (9),  their  first  order  expli¬ 
cit  approximations,  given  by  expressions  similar  to  (8),  will 
be  of  high  quality  only  if  the  structural  members  behave  mainly 
in  extension.  On  the  other  hand,  if  the  bending  behavior  is 
dominant,  it  is  better  to  adopt  the  change  of  variables  (24), 
yielding  first  order  explicit  approximations  of  the  form  (2'i) 
(with  p  =  3).  As  a  matter  of  fact,  the  true  situation  is  usu¬ 
ally  a  combination  of  extension  and  bending.  In  a  practical 
structure,  some  members  work  mainly  in  extension,  some  in  flexion 
anil  others,  simultaneously  in  flexion  and  extension.  Whence 
the  idea  of  using  the  following  explicit  approximations,  which 
should  be  valid  in  any  situation  : 
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where  the  coefficients  efi^  and  c  f  ? 3  are  considered  as  constant 

i  J  t  .1 

throughout  the  redesign  phase. 


Because  it  is  no  longer  possible  to  select  appropriate  in¬ 
termediate  design  variables,  the  explicit  approximations  (33) 
cannot  be  obtained  by  merely  using  first  order  lav  I  or  series 
as  in  the  case  of  pure  bending  elements.  However  an  essential 
requirement  is  that  these  approximations  remain  correct  up  to 


J  . 


Che  first  order,  despite  the  fact  that  they  do  not  result  from 
a  strict  linearization  process.  In  other  words  the  following 
equality  roust  hold  : 
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This  condition  insures  that,  at  the  optimum,  the  solution  to 
the  explicit  approximate  problem  satisfies  the  (first  order) 
optimality  conditions  of  the  real  problem,  that  is,  the  appro¬ 
ximate  and  real  restraint  surfaces  have  the  same  tangent  plane 
(see  Fig.  1).  As  a  result,  the  GOC  approach  should  converge 
to  a  true  (at  least  local)  minimum  weight  design.  It  will  be 
shown  in  this  section  how  such  first  order  explicit  approxima¬ 
tions  can  be  obtained  for  various  types  of  behavior  constraints 
and  structural  models. 


4.1.  Displacement  Constraints 


The  key  idea  to  obtain  explicit  approximations  of  the  dis¬ 
placement  constraints  is  to  come  back  to  the  virtual  load  pro¬ 
cedure,  which  permits  decomposing  any  static  response  quantity 
into  the  contributions  of  each  element.  The  expression  (b)  can 
be  rewritten  in  the  more  general  form 
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Now,  for  a  rather  general  class  of  structural  models,  each  ele¬ 
ment  stiffness  matrix  can  be  assumed  to  have  the  following  ex¬ 
plicit  form  in  terms  of  its  design  variable  |  14  | 
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where  the  matrices  K^P^  are  independent  of  the 

Note  that  most  often,  at  least  one  of  the  K  f 1  ^ 
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zero  ( for  example  K.  is  zero  in  the  stiffnes 

a  flat  shell  element).  1  n  t  r  oduc  i  rip,  (17)  into 

tli  at  a  convenient  explicit  approximation  of  a 

traint  (5)  is  : 
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where  the  coefficients 


c<p)  =  ( qT  K<p)  q.)  a2p 
lj  i  i 


(38) 


(  3  9  ) 


are  assumed  to  be  constant  during 
dient  of  this  explicit  approximate 
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On  the  other  hand,  differentiating  (37)  and  inserting  the  result 
into  (36)  shows  that 


P“l 


P 


T 


q 


K.(P) 


(4  1) 


Therefore  it  can  be  concluded  that  the  expressions  (38)  represent 
first  order  explicit  approximations,  in  that  they  restitute  the 
exact  values  of  the  constraints  and  their  first  partial  deriva¬ 
tives  at  the  design  point  a°  where  the  structural  analysis  is 
made  : 
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An  alternative  approach  is  to  employ  simple  linear  approxi¬ 
mations  in  the  reciprocal  design  variables  (9),  i.e.  the  first 
order  Taylor  series  expansions  (10).  Then  the  explicit  sub  pro¬ 
blem  exhibits  exactly  the  same  form  (14-16)  as  in  the  case  of 
thin  walled  structures  and  dual  methods  need  not  be  modified  as 
they  will  be  in  section  5.  Of  course,  as  previously  mentioned, 
the  convergence  of  the  overall  optimization  process  might  be 
lowered,  or  even  become  unstable,  in  the  case  where  the  bending 
behavior  is  dominant  in  most  of  the  elements.  The  reader  is  re¬ 
ferred  to  sections  6.1  and  6.5  for  numerical  examples  comparing 
the  explicit  approximations  (8),  (25)  and  (38). 

4.2.  Stress  Constraints 


The  situation  is  much  more  delicate  for  stress  constraints 
than  for  disp1 acement  constraints.  In  contrast  to  the  case  of 
thin  walled  structures  modeled  by  bar  and  membrane  elements,  the 
stress  matrices  are  no  longer  constant  in  the  case  of  flexion- 
extension  elements.  For  illustration  consider  again  a  fiat  shell 
element.  The  stresses  in  the  upper  sheet  can  be  computed  in 
terms  of  the  generalized  displacements  by 
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where  o,  ,  are  matrix  representations  ol  the  total 

stress,  membrane  stress  and  bending  stress  tensors  in  member  k, 
and  T ^ are  the  corresponding  stress  matrices,  which 
are  independent  of  the  design  variables.  Each  displacement  com¬ 
ponent  admits  a  first  order  explicit  approximation  of  the  form 
(see  Eq.  33) 
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Therefore,  a  natural  choice  for  the  approximation  of  a  stress 
component  could  be  as  follows  : 
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The  expression  (.43)  is  still  a  first  order  approximation. 

It  is  no  longer  separable,  but  keeps  a  simple  algebraic  form 
suitable  to  specialized  algorithms  such  as  those  based  on  dual 
methods.  However,  when  considering  the  true  stress  limitation 
(e.g.  upper  limit  on  the  Von  Mises  stress),  the  explicit  cons¬ 
traint  becomes  much  more  complicate.  For  these  reasons,  and 
before  finding  something  better,  it  has  been  decided  to  follow 
the  alternative  strategy  proposed  at  the  end  of  section  4.1, 
that  is,  to  employ  simple  linear  approximations  of  the  form  (10). 
To  compute  the  gradient  of  the  stress  constraint  with  respect  to 
the  reciprocal  variables  (9),  we  note  that,  from  the  equivalent 
Von  Mises  stress 
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2  x .  x  y  3  x .  y  x  i  x . 
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I  /  2 

k 
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(46) 
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xy 


xy  3xi  '  k 


1  ,  (4 7  ) 


On  the  other  hand,  from  (44)  it  follows  that  Lin-  derivative  o 

as 

(f  ) 


any  stress  component  o  or  i  lias  the  f  o  mi 


3  o ,  u 

_ x  =  T  qq  _  i 

3  x .  k3x.  ^  i  k  2 

r  1  x . 

l 

where  6^  is  the  Kronecker  symbol. 
4.3.  Frequency  Constraints 


( 4  H  ) 


Constraints  on  natural  frequencies  usually  consist  in  imposing 


1.  .  (  a  ) 


w  .  >  0 

-J 


m 


1  owe  r  limits 


m?  (a) 


J 


(4*0 


f 


1 

I  9  . 


They  are  directly  written  in  terms  of  the  squares  of  the  fre¬ 
quencies,  because  these  quantities  naturally  appear  in  the 
eigenproblem  characterizing  the  structural  modal  analysis  : 


Kq.  -  ut  M  q.  =0 
J  J  J 


(  5(J) 


in  this  equation  K  and  M  represent  the  stiffness  and  mass  matri¬ 
ces,  and  (q.,  j  =  1  , m)  are  the  modal  displace  me nts,  i . v .  ,  the 

J  _  2 

eigenvectors  solution  of  (50),  associated  with  eigenvalues  .,‘1  . 
The  structural  mass  matrix  has  a  linear  form  in  terms  of  the 
design  variables  : 


n  n 

M  =  . Z ,  M .  +  M  =  .  £  ,  a .  M .  +  M 
i=l  l  c  i=l  i  l  e 


(5  1) 


where  and  are  independent  of  the  design  variables.  M ^ 

denotes  the  mass  matrix  of  the  ith  element  when  a.  =  1.  M  re- 

i  c 

presents  the  contribution  of  the  non- structural  masses,  such  as 
equipments,  fuel,  etc...  It  is  well  known  that  the  first  deri¬ 
vatives  of  the  frequencies  with  respect  to  the  design  variables 
are  given  by  |  see  for  example  Ref.  [  2  ],  section  4.2  ]  : 


Du'.  ,  ...  9 K . 

— 1  =  JL  q!(—± 
3a.  m .  j  J  a . 

l  j  J  l 


2  3Mi 
UJ  j  3a. 


)qi 


(52) 


where  m . 

J 


is  the  generalized  mass  of  the  j th  mode  : 


T 

m  .  =  q  .  M  q  . 

J  J  J 


( 5  5 ) 


The  way  to  derive  first  order  explicit  approximations  oi  the 
frequency  constraints  is  less  apparent  than  for  displacement  cons¬ 
traints.  In  this  report,  guided  by  the  work  done  in  Ref  1  2  | 
(section  4.2),  the  following  decomposition  of  the  eigenvalues  in 
terms  of  the  stiffness  and  mass  contributions  of  each  element  will 
be  used  : 


2 
u  . 
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Wj  (  I 


m  . 
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iii 


q!  (K.  -  m2  M .  ) q  . 
J  i  J  i  J 


(54) 


=  m 
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whe  re 


m  . 
J 


(55) 


represents  the  contribution  of  the  non  structural  masses  to  the 


generalized  mass  m  j  (see  Lq.  hi).  By  taking,  account  of  the  ex¬ 
plicit  definitions  (37)  o l  the  stiffness  matrices  K  ^  and  (hi) 
of  the  mass  matrices  M.  ,  the  ii  i  gh  quality  explicit  approximations 
of  the  frequency  constraints  take  the  form  (>!)),  with 


7  m  i 

u.  =  at  (I  +  -J-) 

.1  J  raj 


c<p) 

i  J 


T  M . 

J  t 


lp)qi 


where  =  1  only  If  p  =  1  aid  is  0  otherwise. 


The  coefficients  c 


(P) 


and  the  modified  limits  u.  are 


i  J  J 

their  values  at  the  current  design  point.  lust  as  fo 

placement  constraints,  it  is  easily  verified  that  the 

(38)  remain  first  order  explicit  approximations  satis 

43).  As  a  matter  of  fact,  they  can  be  interpreted  as 

Taylor  series  expansions  in  terms  of  1  /  x  ^  ,  )/x7,  1 / x ! 

as  independent  variables. 


( 5  f>  ) 


1,2,3  07) 

frozen  to 
r  the  d  i  s  - 
express!  mis 
f  y  i  n  g  (42, 
first  order 
considered 


4.4.  Buckling  Constraints 


Just  as  the  natural  frequencies,  the  critical  load  factors 
Xj  are  defined  through  an  eigenproblcm 


K  q  .  -  '  -  S  q  .  -  0 

J  J  J 


(  3  8  ) 


where  S  represents  the  geometric  stiifness  matrix  and  ( q - , 
j  =  I, m)  denote  tiie  eigenvectors  solution  ot  problem  (58),  asso¬ 
ciated  with  eigenvalues  3  ^  .  I  he  pli  v  s  t  c  a  l  meaning  of  the  q  •  1  s 
is  that  of  displacements  in  the  jth  buckling  node  ,  for  a  criti¬ 
cal  load  factor  The  buckling  constraints  consist  in  imposing 

lower  limits  on  the  buckling  loads  : 


i  .  >  >  . 

.1  I 


j  =  1  ,  m 


(  )  <)  ) 


In  this  work,  the  following  form  oi  t h 
adopted  : 


r.ii  nts  will  b  i 


hL 


j 


2  I  . 

;  ,  .  1  1  . 

hj(a)  —  -  —  ( - j  ?  U  j  =  1,  m  (f-'j ) 

because  it.  has  been  found  that  better  explicit  approximations 

are  generated  when  expanding  the  reciprocal  of  the  buckling. 

loads  rather  then  the  A.'s  themselves. 

J 

The  stiffness  matrix  K  has  the  form  (23).  The  geometric 
stiffness  matrix  is  related  to  the  initial  stress  state  in  the 
elements  and  therefore  it  depends  implicitly  on  all  the  design 
variables  : 

n 

S  =  L I ,  S  £ ( a )  (61) 

It  is  worth  recalling  that  the  matrices  S  ^  are  independent  of 
the  design  variables  for  a  statically  determinate  structure. 


As  explained  in  section  5.2  of  Ref  [1 
tives  of  the  buckling  loads  are  given  by 


the  tirst  d  e  r  i  v  a  - 


J  \  . 

_ J.  = 

3  a  . 
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q  .  S  q  . 
J  J 


T  /jKi  _  ,  3  S 

q  j  3a.  j  3a. 


)q  : 


(62) 


In  opposition  with  the  static  and  dynamic  cases  previously  dis¬ 
cussed,  the  derivatives  appearing  in  (62)  are  not  directly  a  - 
vai table,  because  the  elements  of  the  geometrical  stiffness  ma¬ 
trix  are  function  of  the  stresses  acting  in  the  prebukling  state 


Howe.'er,  by  assuming  that  the  terms  - 


as 

3a; 


are  negligible,  the  gra¬ 


dients  (62)  become  easily  computable.  This  assumption,  which  is 
typical  of  optimality  criteria  approaches  for  static  constraints 
amounts  to  not  taking  into  account  the  effects  of  structural 
redundancy  : 


3  S 

3  a . 

i 


=  0 


i  =  1 


(63) 


in  this  re port, the  following  decomposition  of  the  reciprocal 
buckling  loads  in  terms  of  the  contribution  of  each  element  will 
be  used  : 


—  .  .  q  .  K  .  q 

S  .  1=1  1  1  1 

.1  1 


wh  e  re 


s  ■  -  q  .  S  q 

I  .1 


Subs  Ci  tuc  ill)'  the  explicit  del  ini  t  i  eii  i  )  ;  the  s  r  i  M  ness  is.  a  - 

trices  into  (64)  yields  tlu-  high  «|ua  1  i  ty  «.• p  I  i  .  it  approx  uim  - 

tiuns  of  the  buckling  load  constraints  <  )  under  the  imn  (  Jh  } 

with 


,  .  q'{  K<f0  q  .  , 

cfP)=  J.  —  l  :i~''  ,,  .  ,  : 

i  J  2  r 

s  . 

J  J 

Again,  the  coefficients  c  f  l5  ^  and  the  nodi  tied  l  i  n  i  t  s  a.  are 

^  J  } 

frozen  to  their  values  at  the  current  design  point.  lust  e, 
for  the  displacement  and  frequency  constraints,  it  is  easily  vc 
rifieJ  that  the  expressions  (  18  )  are  t  i  r  :>  t  ..  rdc  s  expli.it  a  ppr»- 
xitna  tions  of  the  constraints  (nil  j  ,  s  a  t  i  s  t  •.  i  n  t  • .  ,  a  1 )  ,  provided 
that  the  assumption  (I>1)  can  be  int  reduced  into  the  gradient 
(62)  (statical  de  L  e  rin  i  n  a  cy  )  . 


The  reason  why  the  buckling  n  a  i  a  t  *  ire  treated  bv  e  x  ■ 


panding  tlie  reciprocal  critical  loads 


!  <i  .  6  4  )  is  that,  to  r 


a  thin  walled  structure,  the  a  p  p  r  o  i  n  a  t  i  i  r  t  1  8  )  takes  the  lorm 

(  8  )  ,  where  the  coefficients  c  .  .  r  e  n  a  i  n  n  nsu.it  along  the  s  e  a  1  i  n  , 

1  ! 

line.  This  property,  which  is  well  detailed  in  section  1  o  t 
Ret  |  2  |  ,  lias  been  found  very  important  tor  the  ease  of  stress 

and  displacement  constraints.  (.‘cone  t  r  i  c  a  1  1  v  it  means  that  the 

real  restraint  surlacu  ii  .  ( a )  =  ()  is  re  p  lac  id  bv  a  tangent  snr- 

'v  .  }  . 

face  hj(a)  =  0  at  its  point  of  intersect  ion  with  the  scaling 
line  (see  K i g .  i  of  K e i  .  |  2  |  )  . 


3.  SULUTLON  OF  TIIF.  EXPLICIT  PROBLEM 


From  the  foregoing  developments,  it  appears  t  li  a  t  ,  at  eaeli 
stage  o  t  tile  optimization  process,  the  following.  s  ub  p  r  oh  I  t-m  o'  1 1 : .  t 
be  solved 


minimize 
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subject  to 
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1  ph 

u  . 
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^  a  . 
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(70) 

Unless  the  summation  on  p  in  (69)  is  restricted  to  a  single 
value  of  p  (i.e.  approximation  of  the  form  (25)  with  p  =  1,  2 
or  3),  it  is  no  longer  possibLe  to  find  intermediate  variables 
| i.e.  x^  given  in  (24)  1  in  terms  of  which  the  explicit  cons¬ 
traints  (69)  would  be  linear.  Therefore,  the  primal  solution 
of  problem  (68-70)  is  more  difficult  to  achieve  if  a  gradient 
projection  type  of  algorithm  is  employed  as  in  the  mixed  method 
developed  in  a  previous  work  for  thin  walled  structures  |  1  ,  l  !>  |. 
The  expressions  (69)  are  still  explicit  and  they  continue  to  ex¬ 
hibit  a  simple  algebraic  form.  Consequently  a  general  purpose  op¬ 
timization  a  1  go  r  i  thin  s  uch  as  NKWSHMT  |  1  6  |  could  ensile  be  adapted  to 
take  the  constraints  (69)  into  account.  However,  because  they 
are.  still  additively  separable,  resorting  to  dual  methods  remains 
probably  the  best  strategy,  just  as  in  the  case  ot  thin  walled 
structures  [1,  2,  6,  8  ]. 

The  minimization  problem  (68-70)  can  be  solved  efficiently 
as  an  auxiliary  maximization  problem  in  the  m  1  a  g  r  a n  g i a n  multi¬ 
pliers  r j  associated  with  the  explicit  behavior  constraints  (b9). 
This  dual  problem  reads  as  follow  |  8  ]  : 

n  nt 


maximize  H  (  r )  =  . L .  e .  a.(r)  +  . I ,  r.  g  .  ( r )  (71) 

i  =  l  t  i  j  =  1  J  j 
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de  no  t e 

the  components 

o  f 
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1  i  <■  n  t  , 

w  h  i  c  h 
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primal  constraints  : 
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(  7  i) 


The  primal  variables  .  i  .  (  r)  a  i  ••  related  t  a  t  1 ;  <  •  dual  variables  r  . 

L  1 

through  the  loll  viw  i  ;i-,  ■  j  u  a  s  i  -  u  u  e  ■  n  s  i  r  a  i  in.  .1  mi  ni  mi  /.at  inn  p  r  o  b  I  .  ■  m 
(see  fin  example  Ke!  j  n  j  )  : 


m  l  n  l  ni  r  z  e 


I  ■  1  i  ;  - 


<  p  ) 
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,  P 


u j  )  (  74  ) 


subject  to  a. 


a .  v  a 

l  i 


(  7  i)  ) 


Because  of  the  separability  id  t  his.  problem,  it  can  be  decompose 
into  u  one - d i me ns i ona i  minimization  problems  m  the  form 

3  <i'° 

min  I  <  i  <i  ^  -*  ,  - -  |  (76) 
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where 


(  p )  ‘.n 

C : 1  =  p  ,  r  .  c: 

i  J  =  1  J  i .! 


(P) 


(77) 


Setting  to  zero  the  first  derivative  ol  tie-  single  variable 
function  appearing  in  {lb),  it  is  .ecu  that  the  h  "active"  de¬ 
sign  variables  can  be  obtained  b  solving  the  nonlinear  alge¬ 
braic  e  illations 


:>  ‘i 


(  V  ) 


P  =  1  a  P  ’  I 

i 


I  ,  n 


(78) 


Note  however  that  when  solving  be  side  constraint  ( 7 '* ) 

must  be  taken  into  account  and  -bat,  if  several  values  satisfy 
equation  (78),  (he  one  that  mi  n  i  mi  /.i.  s  ti.o  i  unction  (  7 1> )  must  be 
retained.  Standard  techniques  can  be  employed  to  deal  with  such 
a  simple  one  - d i me  ns i ona  i  problem  (see  t. g.  I  17  I  ). 


In  many  cases,  as  previously  neniieo.td,  at  least  one  ol  the 
terms  in  the  summation  -n  p  in  (  •'  •>  >  a  n.i  (78)  disappears,  and 
it  is  possible  to  treat  (he  problem  analytically.  As  an  illus¬ 
tration,  consider  again  the  ease  o I  a  !  1  a  t  shell  element,  where 
the  term  in  p  -  l  is  missing, 
form  (omitting  the  index  i)  : 


i‘  rub  1  e  m 


(7b)  exhibits  then  the 


I', . 


1 1 ) 
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(  /  9  ' 


while  equation  (ft-)  i.;:  i !  n  I  !i’u' 

.  (  i  i  i  i  ) 

'  <  (HO) 


or,  sotting 


i  1  ■■ 

C  s  0 


(8  1) 

(82) 


The  constant  4  is  known  to  be  i  live,  and  ovt  ry  thing  depends 

thus  on  the  sign  of  the  con  .tints  2.  ^  and  0  ^  ^  ,  which  repre¬ 
sent  the  contribution  ot  the  membrane  and  the  bending,  respec¬ 
tively.  In  general,  equation  (82)  admits  two  solutions,  and  the 
one  that  renders  minimal  the  function  (,/*’)  must  be  adopted. 


From  the  foregoing  developments,  it  appears  that  the  dual 
function  (71)  can  be  c  ■  ■  tin  i  d  e  i  d  as  i  function  oi  the  dual  vari¬ 
ables  only.  The  dual  problem  stated  in  (71,  72)  exhibits  an 

attractive  feature,  name  1 y ,  it  is  a  quas i -  uncons t  rained  problem, 
because  taking  care  ot  fin  no uuo ga t i v i f y  constraints  (72)  on  the 
dual  variables  is  stiai  nti  >  v  •  .  i :  i  .  ■  u  ■■me  iioiuug  a  tivc  dual 

variables,  the  c  o  r  t  e  s  p  o  n  u  ■  n  e  pi  i  .m  i  v .  i)  i  i  r.  J  e  are  computed  by  sol¬ 
ving  (7b)  or  (  /  8  >  ind  the  :>  r  i  m..  j  ;  .  I  i  a  i  u  t  s  a  re  evaluated  by  u- 
s  i  ng  ( 7  J  )  .  The  dual  f'uiie  ».  inn  !  .  i  i  aid  its  gradient  (73)  are  then 
directly  known  and  a  !  e  a.  s  i  b  i  o  a  ..  .•  o  n  t  d  i  recti  on  can  therefore  be 
determined.  In  a  secctul  »rJci  algorithm,  the  hessian  matrix  of  the 


dual  function  has  to  hi-  computed.: 
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and  by  differentiating  the  definition  (78)  ol  a . ( r )  with  respect 

to  r ,  ,  it  follows  that  . 
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On  tlie  other  hand,  i  i  is  oas  1  1  v  see:,  f  ruin  (77)  th.it 


Regrouping  terms,  it  !  in.ii  1)  he  concluded  that  the  terms  of 

the  hessian  matrix  « r  e  given  h  v 
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where  the  summation  or:  t  is  restricted  to  the  n  active  primal 

variables,  that  as,  thogr  that  arc  not  lixed  to  a  lower  or  an 

j  a  • 

upper  limit  (of  course  - - =  0  for  a  passive  variable).  Knowing 

‘ 1  k 

the  gradient  (73)  and  tiie  hessian  matrix  (87)  furnishes  the 
Newton  searcli  direction 


=  -  H 


The  next  dual  point  is  then 


given  hv 


where  a  is  the  step  1  <■  tig  -J:  a 
often  a  regular  Newton  ....  a 
lected,  however,  the  va  in,.  .  t 


a  1  c  a  ;  n  e  d  i  r  e  t  t  i  o  n  /  .  M"  s  t 
:  ■  :  c  r  i  i  in  (80)1  is  sc 

:>i  »  :  I'.' I  t  i  me  s  be  lowered  to 


prevent  one  of  the  dual  variable  ..  i  ram  he  coming  negative  |  6  | 


The  second  order  dual  opt;  mix,  r  i  tap  i  (  mo  n  t  e  d  in  SAMO.KF  |  IK  | 
has  been  especially  devised  so  Lti.ai  ir  seeks  the  maximum  ot  the 
dual  function  by  operating  in  a  s  e  i|  ue  n  i  <■  of  dual  subspaces  w  i  t  h 
gradually  increasing  d  i  me  psiiai.  In  this  way,  the  effective  di¬ 
mensionality  ot  the  maximisation  problem  never  exceeds  the  num¬ 
ber  of  active  Idi.ivi  rr  on s  t  i  1 1  n  Is  ,  which  correspond  to  non¬ 

zero  dual  variables.  Past  experience  with  thin  walled  structures 
indicate  that  this  number  is.  relatively  small  in  practice,  w  h  i  c  li 
explains  the  remarkable  el  i  i  >  i  emv  of  the  dual  method  |  6  j  or 
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con  J  urdt  r  <1  u  i  '  i  t  f  •  •  i  i  t  .1 ,  first  introduced  in 
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6.  NUMERICAL  APPLICATIONS 


All  the  examples  presented  involve  rather  sophisticated 
fiat  shell  elements  that  are  characterized  by  a  displacement 
field  cubic  in  extension  and  quintic  in  flexion  (hybrid  qua¬ 
drangular  flat  shell)  [  18  ]  . 

6.1.  Cantilever  Beam  with  End  Moment 


The  first  example  is  concerned  with  a  cantilever  beam  loaded 
with  a  concentrated  moment  at  its  free  end  1  see  Fig.  2  ] .  This 
problem  has  been  previously  solved  numerically  by  PRASAD  and 
HAFTKA  |  19  ]  using  an  extended  interior  penalty  function  for¬ 
mulation  for  a  beam  having  the  following  properties  :  length  = 

10  in,  width  =  1  in,  applied  moment  =  540  in. lb  ,  Young's  modu¬ 
lus  =  107  psi,  Poisson's  ratio  =  0.3  and  mass  density  =  0.3  lb/in 
There  exists  also  an  analytical  solution,  which  was  obtained  by 
HAUG  |  20  ]  .  The  analytical  optimum  design  for  a  displacement 
limit  of  0.5  in  at  the  free  end  and  specified  allowable  stresses 
of  30,000  psi  is  as  follows  (see  Fig.  2)  : 


a(x)  =  0.30 

for 

0  <•  x  <  2 .  3 

(90) 

a(x)  =  0.244  ^ 

for 

2. 3  <  x  <  10 

(9  1  ) 

It  is  worth  noticing  that  the  beam  being  statically  determinate, 
the  redesign  relation  (31)  is  exact  (with  p  =  3)  and  reduces  to 
(90,  91).  Therefore  the  optimal  solution  should  be  generated 

after  one  structural  analysis  only. 

The  structure  is  discretized  using  20  quadrilateral  plate 
bending  finite  elements  as  indicated  in  Fig.  2.  There  is  no 
design  variable  linking  and  the  problem  involves  thus  20  inde¬ 
pendent  design  variables.  The  problem  was  first  solved  by  using 
the  cubic  expansions  (25)  with  p  =  3  and  the  dual  algorithm  des¬ 
cribed  in  section  5.  As  previously  stated,  because  the  structure 
is  statically  determinate,  the  explicit  problem  (68-70)  is  exact 
and  the  optimum  design  is  obtained  in  one  single  analysis  what¬ 
ever  may  be  the  initial  thickness  a°.  Then,  the  same  example 
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The  second  example  cons i  :  in  minimizing  the  weight  of  the 

simply  supported  square  plate  shown  on  !■  i  g .  3.  It  is  subjected 

to  a  concentrated  load  o  i"  )0U0  its  center,  where  the  de¬ 

flection  is  limited  to  0.02  m .  :,\  svi.-me  try  only  a  quarter  of 
tlu*  plate  has  to  be  analyzed.  >  he  mesh  involves  2  5  plate  ele¬ 
ments  and  175  degrees  of  freedom.  iht*  material  properties  are 
as  follows  :  Young's  modulus  !•  .  .in*'  N  m'  ,  Poisson's  ratio 

v  =  0.3  and  weight  density  '•*.  ■  Vhe  dimensions  of 
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p  =  2  in  the  redesign  relations  (11),  almost  50  iterations  were 
r  eq  u  i  re  d  . 


b . J .  Simply  Suppot  ted  Square  Plate  with  Frequency  Constraint 

Attention  is  now  incused  on  the  simply  supported  square  plate 
shown  in  Fig.  4.  This  example  is  taken  from  Ref  (21  j  and  it  is 


30  . 

concerned  with  the  minimum  weight  design  of  the  plate  subject 
to  a  natural  frequency  constraint.  The  dimensions  of  the  plate 
are  10  in  by  10  in  and  its  material  properties  are  Young's  mo¬ 
dulus  =  30  x  1  0  ^  lb/in^  and  Poisson's  ratio  =  0.3.  In  Ref  I  31  | 
i't  is  stated  that  the  uniform  plate  with  pm  =  1400  is  taken  as 
the  initial  estimate  to  the  optimization  problem,  where  p  is  the 
mass  density  and  w,the  minimum  fundamental  frequency.  Therefore, 

in  the  present  study,  assuming  steel  material  with  miss  density 
3 

p  =  0.283  lb/in  ,  the  minimum  frequency  was  chosen  as  =  u>/2  n  - 
11.2  Hz.  Iteration  history  data  are  presented  in  Table  4  for 
three  cases  differing  by  the  minimum  thickness  constraint  (a  =  0. 1 
0.05  and  0.001  in).  The  initial  design  in  cases  1  and  3  corres¬ 
ponds  to  a  uniform  thickness  a°  =  0.2  in,  while  in  case  2,  n° 
equals  0.12  in.  Final  uesignsare  illustrated  in  Fig.  4.  By 
symmetry  only  a  quadrant  of  the  plate  has  to  be  analyzed  and 
designed.  It  can  be  observed  that  the  design  obtained  in  Kef 
[21  ]  is  different  from  the  design  achieved  in  this  study  (case  I) 
However  both  designs  have  about  the  same  weight  (0.765  lb  and 
0.752  lb).  It  is  also  worth  pointing  out  that  the  design  of  pla¬ 
tes  in  bending  is  an  analytically  difficult  problem,  leading  to 
multiple  local  optima  and  suggesting  that  the  optimal  plate  should 
be  made  up  of  an  infinite  number  of  stiffeners  |  22  )  . 

6.4.  I-beam  Structure 

The  third  example  involves  the  I-beam  structure  schematized 
in  Fig.  5.  The  problem  consists  in  minimizing  the  weight  ol  the 
beam  while  imposing  lower  bounds  on  the  frequencies  of  the  three 
first  eigenmodes  :  flange  flexion,  torsion  and  web  flexion.  De¬ 
tailed  data  can  be  found  in  Ref.  [  2  ] .  In  a  first  optimization 
exercise,  a  pure  membrane  model  was  employed.  It  involves  35 
second  degree  displacement  elements,  including  10  fictitious 
diaphragms  (without  masses).  These  dummy  members  are  introduced 
to  obtain  a  satisfactory  representation  of  the  torsional  mode. 

Only  5  analyses  are  sufficient  to  generate  an  optimum  design  for 
this  membrane  model. 


3  I  . 

However,  when  this  final  design  was  analyzed  by  using  a  more 
accurate  model  made  up  of  flat  shell  elements,  the  torsional  fre¬ 
quency  (inode  2)  was  seen  to  be  violated  by  10  7  .  Therefore  tin 
problem  was  again  solved  with  this  u  e  w  m  o  d  <  I  ,  by  • <  s  o  r  t  i  m*  ' 1 '  the 
theory  proposed  in  section  4.3  and  to  the  dual  optimizer  descri¬ 
bed  in  section  5.  Iteration  history  data  are  illustrated  in 
Fig.  5  and  the  final  designs  are  given  in  Table  3  for  both  finite 
element  models  of  the  I-beam.  It  can  be  seen  that  the  use  of 
flat  shell  elements,  although  yielding  slower  convergence,  gives 
satisfactory  results. 

6.5.  U-beam  Structure 

In  an  attempt  to  consider  a  case  where  both  flexion  and  ex¬ 
tension  loadings  play  an  equally  important  role,  the  U-beam 
structure  depicted  in  Fig.  6  was  optimized.  The  2000  kg  load 
acting  at  the  tipproduces  torsion  of  the  beam.  As  a  result  the 
upper  flange  behaves  mainly  in  flexion,  the  web,  in  extension, 
and  the  lower  flange,  both  in  flexion  and  extension.  For  sim¬ 
plicity,  only  displacement  constraints  are  considered  :  the  tip 
deflection  is  limited  to  ().!  m,  while  the  relative  lateral  dis¬ 
placement  is  limited  to  0.005  m  (see  Fig.  6).  Three  different 
methods  are  employed,  which  ditfei  by  the  explicit  approximations 
used  for  tlie  displacement  constraints  : 

cast;  I  :  full  expansion  (  ,1  '■  )  .  yielding  problem  (68-70  )  with 
p  =  1  and  3  only  (the  lerir  in  p  =  2  being  zero) 

case  2  :  linear  expansion  (.  2  >  witii  p  -  !),  yielding  problem 

(14-16) 

case  3  :  cubic  expansion  (25  with  p  =  3 ) . 

lhc  iteration  history  data  are  given  in  fable  6  and  the  final 
designs  in  Table  7. 

Lt  is  interesting  to  notice  that  the  three  methods  give 
similarly  good  results,  which  suggests  that  the  use  of  simple 


linear  Taylor  series  expansion  in  terms  of  the  reciprocal  design 
variables  might  be  a  good  strategy  for  most  of  the  structural 
optimization  problems.  This  idea  was  previously  stated  by  AUSTIN 
[  23  ]  ,  but  not  believed  valid  by  many  people  when  he  published 
his  paper  ! 


7.  CONCH'S  IONS  AND  Kfci  COMMAN  D  AT  1  UNS 


i  i  . 


It  is  now  widely  recognized  lint  a  powerful  and  rather  ge¬ 
neral  approach  to  structural  optimization  is  achieved  bv  repla¬ 
cing  the  original  problem  with  a  sequence  of  explicit  approxi¬ 
mate  problems.  This  approach  was  initially  conceived  for  thin 
walled  structures  modelled  bv  bar  and  membrane  finite  elements, 
as  well  in  the  context  of  optimality  criteria  techniques  as  in 
the  framework  of  mathematical  programming  methods  using  approxi¬ 
mation  concepts.  It  has  been  extended  in  this  work  to  deal  with 
structural  systems  made  up  of  beam,  plate  and  flat  shell  elements, 
with  behavior  constraints  placed  on  displacements,  stresses,  na¬ 
tural  frequencies  and  critical  buckling  loads.  The  method  p re¬ 
sented  uses  a  second  order  dual  algorithm  to  solve  each  explicit 
sub  problem.  The  convergence  properties  are  independent  of  the 
number  of  design  variables,  which  is  typical  of  optimality  cri¬ 
teria  types  of  approach  as  well  as  linearization  techniques  in 
mathematical  programming.  As  a  result  large  structural  systems 
can  be  treated  at  the  expense  of  a  few  finite  element  analyst's. 

Several  difficult  points  remain  to  be  clarified.  First  it 
is  not  well  known  whether  the  explicit  problem  can  be  solved 
conveniently  in  any  case,  because  its  lack  of  convexity  might 
lead  to  discontinuity  in  the  dual  function  gradient.  Secondly, 
the  use  of  simpLe  linear  expansions  with  respect  to  the  recipro¬ 
cals  of  the  element  transverse  sizes,  gives  rise  to  satisfactory 
results  in  many  cases.  This  suggests  that  more  complicate  expli¬ 
cit  approximations  such  as  those  proposed  in  this  work  mi  g  lit  not 
be  necessary.  Finally,  there  remains  the  question  of  the  stress 
constraints,  for  which  it  is  difficult  to  generate  adequate  first 
order  explicit  approximations.  Also  s t r e s s - r a t i o i ng  algorithms 
should  be  considered  lor  extension- flexion  elements,  because  they 
correspond  to  much  more  inexpensive  zero  order  approximations. 
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Table  2  Final  Designs  for  Cantilever  Beam 


E  lenient 
No. 


Analytical 
Ref  [  20] 


Thick  ness  (in) 


Ref  [  19  ] 


Linear 
Expans  ion 


ub  i  c 
Expansion 


1 

2 

3 

4 

5 

6 

7 

8 
9 

1  0 
1  1 
1  2 
1  3 
14 
1  3 
1  6 
1  7 
1  8 
1  9 
20 


Weight 

(lb) 


0 . 3 
0.  3 
0. 3 
0 . 3 
0 . 3 

0.3142 
0.3276 
0. 3396 
0.3503 
0.3602 
0.3693 
0.3778 
0 .3858 
0.3933 
0 . 4004 
0.407  1 
0.4135 
0.4197 
0.4255 
0.43  12 


1  .0750 


0 . 3004 
0 . 3006 
0.3018 
0 . 3056 
0.3109 
0.3167 
0.3231 
0 . 3298 
0 .3367 
0.3437 
0. 3508 
0.3579 
0 .3652 
0 .3726 
0 . 3803 
0.3881 
0 . 3960 
0.4040 
0.4119 
0.4197 


1  .0453 


0.  3 
0.  3 
0.  3 
0 . 3 
0.  3 

0.3167 
0.3287 
0.3327 
0.3350 
0.3379 
0. 3450 
0.3574 
0.3711 
0.3851 
0. 3989 
0.4121 
0.4246 
0.4358 
0.4446 
0.4533 


1.0768 


0 . 3 
0  .  3 
0 . 3 
0  .  3 
0  .  3 

0.3136 
0.3270 
0. 3389 
0.3497 
0.3596 
0.3687 
0.3771 
0.3851 
0.3926 
0.3996 
0.4063 
0.4126 
0.4182 
0.4225 
0 .4268 


1.0797 


No  •  of 
Ana  lyses 


/ 


8 


2 


Table  4  : 


Iteration  History  Data  foi  Square  Plate 
with  Frequency  Constraint 


Ana  lysis 

No. 

Case  1 

‘ 

a  =  0  .  1  in 

Case  2  : 

a  =  0.05  in 

Case  3  : 

a  =  0.001 

Weight 

(lb) 

Frequency 

(Hz) 

Weight 

(lb) 

Frequency 

(Hz) 

Weight 

(lb) 

F  r  e  q  u  o  n  i 
(Hz) 

1 

1.415 

19.6 

0.849 

11.8 

1.415 

i  y .  5 

2 

1.097 

16.0 

0 . 786 

11.5 

1  .097 

16.0 

3 

0.913 

13.6 

0. 750 

11.5 

0 . 902 

13.7 

4 

0.812 

12.2 

0. 708 

11.7 

0.778 

12.4 

5 

0.771 

11.5 

0.658 

11.6 

0.696 

11.9 

b 

0.757 

11.3 

0.63  1 

11.4 

0.647 

11.5 

7 

0 . 753 

11.2 

0.617 

11.3 

0.622 

11.4 

8 

0.752 

1  1  2 

0.605 

11.3 

0 . 604 

11.4 

9 

0.752 

11.2 

0.596 

11.3 

0.588 

11.4 

1  0 

0.589 

11.3 

0.5  74 

11.3 

1  1 

0.583 

11.3 

0.559 

11.4 

1  2 

0.577 

11.3 

0.542 

11.4 

1  3 

0.52  1 

11.3 

IBM  370-158 

CPU  time 

2  1  8 

276 

307 

(sec) 

_ 

nge  d 


Table  5  :  Final  Designs  fur  I-beam  Structure 


thickness  (mm) 


membrane  model 
flat  shell  model 


Upper  Flange 


15.80 

17.30 

11.67 

6.143 

1.813 

16.59 

15.09 

10.08 

5  .  090 

1  .36  3 

Web 


5.101 

3 .602 

3.329 

3.294 

1  .  9  9  7 

2 .056 

4.173 

3.965 

3.936 

2.313 

f 

S  up ported 


Lower  Flange 

I - 1 - 1 - - - 1 


1.79 


15.40 
20.5  1 


17.01 

17.98 


11.55 

12.21 


6.074 

6.480 


1.61  () 


Table  6  :  Iteration  History  Data  for  U-Beam  Structure 


1  COMPARISON  OF  REAL  AND 
no.  \  APPROXIMATE  restraint  surfaces 


finite  element  model 

;  i" 

«* -  10"-  -  * 

20  bending  plate  elements 
220  degrees  of  freedom 

FIG.  2  CANTILEVER  BEAM  WITH  END  MOMENT 


final  design  (1/4)  thicknesses  (m) 


symmetrical 


SIMPLY  SUPPORTED  SQUARE  PLATE 
WITH  DEFLECTION  CONSTRAINT 


